CLAIMS 

1. A method for automatically detecting and tracking the 
contour of a starting image, said starting image being 
formed by an array of pixels and each pixel having an 
intensity, being n and m coordinates of a generic pixel, 
comprising the steps of: 

filtering said starting image through an absolute 
central moment e(n,m) of the intensity of the pixels of said 
starting image, wherein said absolute central moment is 
obtained with the following steps: 

- determining for each n,m a local mean calculated on 
a neighborhood about a pixel of coordinates n,m of the 
starting image, obtaining a first filtered image; 

- determining for each n,m a sum of absolute 
differences between the intensity of the pixel having 
coordinates n,m of the first filtered image and the 
intensity of all the pixels contained in a neighborhood 
about said pixel of coordinates n,m of either said 
starting image or a second filtered image derived from 
said starting image, 

- wherein said sum of absolute differences is split 
calculating a sum of positive differences, or positive 
deviation, and a sum of negative differences, or 
negative deviation. 

2. Method according to claim 1, wherein said sum of 
absolute differences is calculated computing the 
differences between said first filtered image and said 
second filtered image, wherein said second filtered image 
is obtained for each n,m from a local mean calculated on a 
neighborhood about the pixel of coordinates n,m of said 
starting image. 

3. Method according to claim 1, wherein said absolute 
central moment e(n,m) is calculated in a generalized way as 
follows : 
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where : 

n,m are the coordinates of a pixel of a map f(n,m) of 
said image; 

w^(n,m), W2(n,m), W2(n,m) and w^(n,m) are four weight 
functions defined on four circular domains 0i, ©21 
©3 and ©4, respectively of radius r^, r^, and 
and defined as ©,. = {(A:,/) g : V^^T? < ; 

® is a convolution operator; 

//j(«,m) = ^ ^f{n,m)w^{n^m) is a mean value on domain ©1 of 

said map and of said first filtered image; 
//2(n5/w) = ^ ^f{n-k^m-l)w2{k^l) the mean value on domain ©2 

of said map and of said second filtered image. 

4. Method according to claims 3 and 4, wherein, starting 
from said absolute generalized central moment, said 
positive deviation ep(n,m) and said negative deviation 
e^(n,m) are used as further filters, defined as 

ep in, m) = (/i, m) ® ^ ^ (n, m)- ju^in-k^m- l))w^ (k, I) 
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(7) 

e„ (w, m) = W4 (w, m) ® ]^ ^ («, m)- ju^in-k.m- l))w^ (k, I) 
where domains and ©3^ are defined as: 

(8) 

©3„ ={(/:,/)€ ©3 \n,{n,m)<^^{n-k,m-l)} 

5. Method according to claim 4, wherein said step of 
computing the absolute generalized central moment of the 
intensity of a pixel comprises the steps of : 

- defining said circular domains ©i, ©2, ©3 and ©4, in a 
neighborhood about each point of the starting image, 
wherein ©1, ©3 and ©4 are centered on n,m and ©2 is a 
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domain centered respectively on each point inside ©3; 

- associating to each domain with i comprised between 
1 and 4, a weight function and computing a mean value 
of a grey levels map for domains ©i and ©2 on the basis of 

5 said w^; 

- computing the absolute generalized central moment e(n,m) 
on the basis of a weight function on domain ©3/ 

- splitting the absolute generalized central moment e(n,m) 
into a positive deviation e^in.m) and a negative deviation 

10 e^(n,m} , thus creating at said discontinuity two partially 
overlapping bell -shaped profiles; 

convoluting the two positive ep(n,m) and negative e^{n,m) 
deviations with weight function W4 on domain ©4; 

- measuring discontinuity using said positive deviation 
15 e^{n,m) and negative deviation e^{n,m) as filters. 

6. Method according to claim 6, wherein said step of 
measuring discontinuity is carried out by tracking of a 
function 

Min (Bp (n,m) , \ e„ {n,m) \ ) 

20 7. Method according to claim 6, wherein said step of 
measuring discontinuity is carried out through a 
subtraction defined as: 
e(n, m)^ep (n, m) ~ e„ («, m) 

8. Method according to claim 6, wherein a DoG filter 
25 (difference of Gaussian curves) is obtained using a sum of 

said positive e^{n,m) and negative e^{n,m) deviations of the 
central absolute moment. 

9. Method according to claim 1, wherein said starting 
images are selected from the group: 

30 - biomedical images, obtained with ultrasonic pulses, PET, 
SPECT, CAT, MR, etc, among which anatomical images, or 
images of function, obtained by means of time sequences of 
anatomical views of a particular zone of an organ, or 
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perfusion images, obtained on the same organ after treatment 
of the patient with substances that enhance the perfusion in 
the organ; images of graphs acquired by a scanner in order 
to convert paper graphs into digital signals. 

5 10. Method according to claim 9 wherein, in case of images 
of graphs, scanned images are filtered with said absolute 
central moment e(n,m) tracking it as a bell -shaped profile 
whose peak is the sought digital signal, a further step 
being provided of computing the digital signal with a 
10 local maxima detection algorithm of standard type. 

11. Method for contour tracking, according to claim 3, 
characterized in that at a discontinuity said generalized 
absolute central moment, calculated as 

e'(«,/w) = ^ Y,\Mi(n,rn)-f(n-k,m'-l)\w,(kJ) (12) 
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15 is compared with a threshold value derived from said 

generalized absolute central moment, calculated as 
' («, m) = w, («, m) ® ^ 2; m) ^f(n-k,m- 1)\ w, (kj) (13) 
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12. Apparatus for contour tracking in video images arranged 
as succession of photograms, characterized in that an 

20 arithmetic logic unit is used and one or more filters are 
used that calculate a positive deviation e^{n,m) and a 
negative deviation e^{n,m) of an absolute generalized 
central moment as defined in claim 1. 

13. Apparatus according to claim 12, wherein said positive 
25 deviation e^{n,m) and said negative deviation e„(/i,/w) are 

defined as: 

(7) 

e„ («, m) = (w, m) ® J] J] U m)- ^^{n-k^m- l))w, {k, I) 
where domains and 03„ are defined as: 
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03p = {(^. /) G 03 : /i, («, m)> n^{n-k,m- /)} 

(8) 

®3n ={ikj)e@,:ju,{n,m)<ju^(n-k,m-l)} 

14. Apparatus according to claim 12 wherein said filters 
consist of four bidimensional convolutors and an 
integrator. 

5 15. Apparatus, according to claim 13, wherein each of said 
bidimensional convolutors is replaced by a cascade of two 
monodimensional convolutors . 



